Draft version May 8, 2009 

Preprint typeset using 1^1^^ style cmulatcapj v. 03/07/07 



THE 2006-2007 ACTIVE PHASE OF ANOMALOUS X-RAY PULSAR 4U 0142+61: RADIATIVE AND TIMING 

CHANGES, BURSTS, AND BURST SPECTRAL FEATURES 

FoTis p. Gavriil^ ^, Rim Dib'', Victoria M. Kaspi'' 

Draft version May 8, 2009 

ABSTRACT 

After at least 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142-1-61 entered an active 
phase in 2006 March that lasted several months and included six X-ray bursts as well as many 
changes in the persistent X-ray emission. The bursts, the first seen from this AXP in >11 years 
of Rossi X-ray Timing Explorer monitoring, all occurred in the interval between 2006 April 6 and 
2007 February 7. The burst durations ranged from 8—3x10'^ s. The first five burst spectra are 
well modeled by blackbodies, with temperatures kT ~ 2 — 6 keV. However, the sixth burst had a 
complicated spectrum that is well characterized by a blackbody plus three emission features whose 
amplitude varied throughout the burst. The most prominent feature was at 14.0 keV. Upon entry 
into the active phase the pulsar showed a significant change in pulse morphology and a likely timing 
glitch. The glitch had a total frequency jump of (1.9±0.4) x 10~^ Hz, which recovered with a decay 
time of 17±2 days by more than the initial jump, implying a net spin-down of the pulsar. We discuss 
these events in the context of the magnetar model. 

Subject headings: — stars: neutron — X-rays: stars — X-rays: bursts — pulsars: individual 
(4U 0142-F61) 



1. INTRODUCTION 

It is now generally accepted that the class of objects 
referred to as "Anomalous X-ray Pulsars" (AXPs) are 
magnetars - young isolate d neutron stars powered by 
, their high magnetic fields (jThompson fc DuncanI 119951 : 

■ [Thompson fc DuncanI [l996l ). AXPs are X-ray pulsars 
with periods in the range 2-12 s, and X-ray luminosities 
(~ lO^'^ — 10'^^ erg s~^) that cannot be accounted for by 
their available spin-down energy. The magnetar model 
was first proposed to explain the dramatic behavior ex- 
hibited by an apparently different object class - the Soft 
Gamma Repeaters (SGRs). SGRs show persistent prop- 
erties similar to AXPs, but they were first discovered by 
their enormous bursts of soft gamma rays (> 10*^ erg) 
and their much more frequent, shorter, and thus less 

■ energetic bursts of hard X-rays. To date, SGR-like X- 
ray bursts have been observed from five AXPs, thus so- 
lidifying t he conne ction between the two source classes 

; (iGavriil et al. 2002: IKasDi et al.llMillWoods et al.ll2005l : 
, iKrimm et~ al. 2006) . For reviews of magn e tar candidate s 
and AXPs sec Woo ds fc ThompsonI (|2006[ ) , iKaspil (|2007D , 
and lMereghettii t\26o^ . 

Thus far, only the magnetar model can ex- 
plain the bursts observed from SGRs and AXPs 
(jThompson fc DuncanI 119951 ). The internal magnetic 
field exerts stresses on the crust which can lead to large 
scale rearrangements of the external field, which we ob- 
serve as giant flares. If the stress is more localized, 
then it can fracture the crust and displace the foot- 
points of the external magnetic field which results in 
short X-ray bursts. The highly twisted internal magnetic 
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field also slowly twists up the external field; magneto- 
spheres of magnetars m ay therefore be globally twisted 
(jThompson et al.]l2002t l. Reconnection in this twisted 
magnetosphere has also been proposed as an a dditional 
mechanism for the short bursts (Lvutikovl l2002f ). 

In addition to bursts, AXPs and SGRs exhibit pulsed 
and persistent flux variations on multiple different 
timescales. An hours-long increase in the pulsed flux 
has been s een to foUow a burst in AXP IE 1048.1-5937 
(IGavriil et al. 2006). On longer timescales, AXPs can 
exhibit abrupt increases in flux which decay on several- 
week timescales. These occur in conjunction with bursts 
and have been suggested as being due to thermal ra- 
diation from the stellar surface after the deposition of 
heat from bursts. Such flux enhancements have been 
observed in SGRs (see IWoods et af] l2001l for exam- 
ple). The flux enhancement of AXP IE 2259-^586 dur- 
ing its 2 002 outburst can also interpreted as burst af- 
terglow (jWoods et al.|[200^ . however, a magnetospheric 
interpretation is also plausible (Zhu et al. 2008). AXP 
IE 1048.1-5937 exhibited three unusual flux 'flares.' In 
the first two, the pulsed flux rose on week-long timescales 
and subsequently decayed back on time scales of months 
(jCavriil fc Kaspil l2004 iTam et al.l l2008f ). These vari- 
ations have been tentatively attributed to twists im- 
planted in the external magnetosphere from stresses 
on the crust imposed by the internal magnetic field. 
AXPs XTE J1810-197, IE 1547.0-5408, and the AXP 
candidate AX J1845-0258 have also e xhibited large 
flux variations [I brahim et al. 2004; Halp ern et al.ll2008l : 
iGottheff fc Vasishtlll998l : lTam et aLl l2006l. however it is 
not clear whether these were of the abrupt rise type as 
in IE 2259-^586 or the slow-rise type as in IE 1048.1- 
5937 CXO J164710. 2-455216 showed a clear abrupt 
rise (jMuno et al.ll2007t ). IRXS J1708-4009 has been ar- 
gued to have flux va riations associated with timing events 
(jlsrael et al.ll2007bf ). Finally, AXP 4U 0142+61 has ex- 
hibited the longest timescale flux variations thus far, in 
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Fig. 1.— Epochs of the RXTE observations of 4U 0142+61 
analyzed in this paper (see § 13.21 1 . Each observation is represented 
by a vertical line. The dashed line marks the entry of the source 
into the active phase. An "x" marks the first observation where the 
pulse profile was significantly different from the long-term average, 
"bl" , "b2" , and "b3" mark the observations containing bursts. To 
the left of the dotted line, the letters refer to groups of observations 
having similar total integration times. To the right, the letters refer 
to groups of observations having similar pulse profiles. 

which the p ulsed flux increa s ed by 29±8% oyer a period 
of 2.6 years (iDib et al.ll2nn7l : iGonzalez et al.llMiflh . 

4U 0142+61 is an 8.7-s AXP. It has a period derivative 
of P = 0.2 X 10^^^, implying a surface dipole magnetic 
field of 1.3 X 10" G. 4U 0142+61 w as monitored by 
RXTE in 1997 and from 2000 to 2007. iGavriil fc Kaspil 
(|2002f) showed t hat 4U 0142+61 generally r o tates with 
high stability. 'Morii, Kawai, fc Shibazakil ()2005f ) re- 
ported a possible timing glitch in 1999 on the basis of 
an Advanced Satellite for Cosmology and Astrophysics 
{ASCA) observation in which the value of the frequency 
w as marginally discrepant \ yith the frequenc y as reported 
by IGavriil fc Kaspil (|2002D . iDib et al.l (|2007D showed that 
the glitch ma y have occurred b ut is not required by the 
existing data. iDib et ahl ()2007l ) also reported on the evo- 
lution of the properties of 4U 0142+61 from 2000 March 
to 2006 April. In particular they reported stable timing, 
and an evolution of the pulse profile in 2—4 keV where 
the dip between the two peak was rising between 2000 
and 2006. They also reported a 29+8% increase in the 
pulsed fiux between 2002 May and 2004 December. As of 
2006 March, in the published fiux history of this source, 
there had been no reports of any X-ray activity such as 
bursts or fiares, as described above for other AXPs. 

Here we report on the first detection of bursts from 
AXP 4U 0142+61, making this the sixth AXP for which 
this phenomenon has been observed. We also report that 
the source appears to have entered an active phase in 
2006 March in which almost every aspect of the emission 
changed. Our observations are described in Section O 
Our burst, pulsed morphology, pulse phase, pulsed flux, 
and timing analysis are presented, respectively, in Sec- 
tions O 112113311331 and EH In Section [H we discuss 
the possible origins of this behavior and the implications 
for the magnetar model. 

2. OBSERVATIONS 

4U 0142+61 has been monitored with Rossi X-ray 
Timing Explorer (RXTE) in 1997 and from 2000 to 2007. 
It has been monitored bi-monthly since 2005 March, with 
a typical observation length of ~5 ks. On 2006 March 
23, the source entered an active phase where many pulsed 
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One PCU switched ON and another one OFF partway through 
the observation, but the total number of active PCUs stayed con- 
stant throughout the observation. 

fiux, spectral, and pulse profile changes were observed. 
We detected 6 bursts in three observations after the en- 
try into the active phase. After each burst detection, 
several RXTE target-of-opportunity (ToO) observations 
were made in addition to the regular monitoring. 

Here, we present a detailed analysis of the three obser- 
vations containing bursts (see Table [T]) . We also present 
an analysis of 91 other observations spanning the 2005 
March 21 to 2007 May 15 time period (MJD 53481 to 
MJD 54235, Observation IDs 91070-05-04-00 to 92006- 
05-32-00). 32 of these observations were prior to the 
entry into the active phase, and the remaining 59 ob- 
servations were after. This long-term analysis was per- 
formed in order to track the evolving pulsed flux and 
timing properties of the source. 

Figure [1] shows a timeline of the 94 analyzed observa- 
tions. Note that two segments of any observation that 
was split counted as separate observations if the seg- 
ments were given different observation IDs. The ranges 
of epochs with an increased density of observations con- 
tain the ToO observations. Prior to the active phase, 
the groups of observations referred to by capital letters 
in the Figure have similar total integration times. The 
groups of observations in the active phase have similar 
pulse profiles. These groups will be referred to in Sec- 
tions 13.21 and 13.51 



3. ANALYSIS AND RESULTS 

All data presented here are from th e Proportional 
Counter Array fPCA: lJahoda et al.ll2006f ) aboard RXTE. 
The PCA is made up of five identical and inde- 
pendent proportional counter units (PCUs). Each 
PCU is a Xenon/methane proportional counter with 
a propane veto layer. The data were collected in ei- 
ther GoodXenonwithPropane or GoodXenon mode which 
record photon arrival times with ~l-/is resolution and 
bins them with 256 spectral channels in the ^^2-60 keV 
band. 

3.1. Burst Analysis 

For each monitoring observation of 4U 0142+61, using 
software that can handle the raw telemetry data, we gen- 
erated 31.25-ms lightcurves using all Xenon layers and 
events in the 2-20 keV band. These lightcurves were 
searched for bursts us ing our automated b urst search al- 
gorithm i ntroduced inlGavriil et al.l (|2002t l and discussed 
further in lGavriil et al.l ( 20041 ) . In an observation on 2006 
April 6, we detected a significant burst, and four more 
bursts were detected in a single observation on 2006 June 
25. The sixth and most energetic burst was detected on 
2007 February 7. The bursts were significant in each ac- 
tive PCU. See Table [1] for the number of active PCUs in 
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each burst observation, as well as for the burst observa- 
tion epochs. 

To analyze these bursts we created event lists in FITS^ 
format using the standard FTOOLS^. For consistency with 
previous analyses of SGR/AXP bursts we extracted 
events in the 2-60 keV band. These events were re- 
duced to the solar sy s tem b arycenter using the position 
found by iPatel et al.l ()2003D for the source. The burst 
hghtcurves are displayed in Figure [2] (LEFT). Our tech- 
niques for characterizing the temporal and spectral prop- 
erties of bursts were discussed in detail in iGavriil et al.l 
1^04) but we repeat them here as some required modi- 
fication in order address the distinct properties of these 
bursts. 

Before measuring any burst parameters we deter- 
mined the instrumental background using the FTDDL 
pcabackest. We extracted a background model 
lightcurve using the appropriate energy band and num- 
ber of PCUs. pcabackest determines the background 
count rate only on 16-s time intervals, so we interpolated 
these values by fitting a polynomial of order 6 to the en- 
tire observation, which yielded a good fit for each data 
set. 

3.1.1. Burst Temporal and Energetic Properties 

The burst peak time {tp), rise time (tr) and peak fiux 
(fv) were de t ermin ed using the methods described in 
IGavriil et all (|2004f ). Usually, to measure the fluence 
for SGR and AXP bursts, we subtract the instrumental 
background for the lightcurve, integrate the light curve 
and fit the result to a step function with a linear term 
whose slope is the "local" background rate. The burst 
fluence in this case is the height of the step function. 
Although this technique worked well for the first burst, 
which was a short, isolated event, it was not appropriate 
for bursts 2,3, and 4 because they had overlapping tails, 
and bursts 5 and 6 had tails that extended beyond the 
end of the observation. Thus, we opted to fit the burst 
flux time series, f(t), with exponential rises and decays: 



fpe^ ) + b{t) + B 



t < t. 



p 



fit) — \ ^(tihL\ _(1z1p.\ 

ifp - ft)e ^ ^ +fte y ) t> tp. 
+bit) + B 

The burst peak time {tp), the burst rise time {tr), the 
peak flux (/p), the instrumental background model {b{t)) 
determined using pcabackest and the "local" back- 
ground {B) were held fixed at the values determined by 
the methods outlined above. We fit for the flux of the 
long-term tail component of the burst, and for the short 
(r/) and long-term (rt) decay times of the tail of the 
burst. We then integrated our events to obtain fluence 
lightcurves. We fit the integral of Eq.[l]to our integrated 
events. The fits were performed over the entire span of 
the observations. Figure [2] (LEFT) shows our observed 
and model burst lightcurves; Figure [2] (RIGHT) shows 
our observed and model burst fluence time series. 

As is done for 7-ray bursts and SGR and AXP bursts, 
we characterized the burst duration by Tgo, the time from 
when 5% to 95% of the total burst counts have been 
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collected. To determine the Tgo duration we integrated 
our burst model and numerically determined the 5% and 
95% time crossings. All burst temporal parameters are 
presented in Table [2l 

3.1.2. Burst Spectral Properties 

Burst spectra were extracted using all the counts 
within their Tgo interval. Background intervals were 
extracted from long, hand-selected, featureless intervals 
prior to the bursts. Response matrices were created us- 
ing the FTOOL pcarsp. The burst spectra were grouped 
such that there were at least 15 counts per bin after back- 
ground subtraction. Burst spectra, background spectra, 
and response matrices were then read into the spectral 
fltting package XSPEC^ vl2.3.1. The spectra were fit to 
photoelectrically absorbed blackbodies usi n g the column 
density found by iDurant fc van KerkwiikI (|2006l ) . Only 
bins in the 2-30 keV band were included in the fits. The 
blackbody model provided an adequate fit for bursts 1 
through 5; model parameters are presented in Tabled 

Extracted over its very long Tgg interval (see Table [2]) , 
the spectrum of Burst 6 was well fit by a simple black- 
body, however, there were hints of spectral features. On 
much shorter timescales, the spectrum of burst 6 was 
not well modeled by any simple continuum model and 
the presence of spectral features is clear (see Fig. [3]). To 
study the temporal evolution of these features we ex- 
tracted 5-s long spectra in steps of 1 s from the peak 
of the burst, and repeated the spectral fitting procedure 
outlined above. Fig. [5] displays the spectral evolution 
of the features. The middle panel is a surface plot of 
the change in after subtracting the continuum model 
(blackbody). A vertical slice in the middle panel corre- 
sponds to a spectrum of the burst extracted over a 5-s 
interval, and these spectra were extracted every 1 s . The 
features showed clear temporal variability and were most 
prominent near the onset of the burst (see Fig. S]). 

We looked closely at the first 5 s of the burst (see 
Fig. [3]), and determined that a simple continuum model 
with less than three spectral features was not a statis- 
tically good representation of the data. A simple pho- 
toelectrically absorbed blackbody yielded a reduced 
of 2.54 for 41 degrees of freedom (DoF). Adding one or 
two Gaussian emission lines improved the fit but not suf- 
ficiently, yielding reduced x^s of 1.89 for 38 DoF, and 
1.67 for 35 DoF, respectively. Only a photoelectrically 
absorbed blackbody plus three Gaussian emission lines 
provided an adequate fit, with reduced 'x^—1.21 for 32 
DoF. With this best-fit model, we find lines at energies 
8.8l|J;2, 14.0t^;3 and 21.8j:? j keV. Our results from fit- 
ting different models to the first 5 s of the burst spec- 
trum are summarized in Figure [3l and the parameters 
returned from our best-fit model are listed in Table [3l 
We note that a power-law with index r=0.9to'2 together 
with three Gaussian emission lines also yielded a good 
fit, with reduced x^=l-20 for 32 DoF. In addition, us- 
ing a power-law instead of a blackbody as the continuum 
model returned values for the Gaussian emission line pa- 
rameters that were consistent with those listed in Table [3] 
at the l-cr level. 

From the surface plot in Fig. 3] (middle panel), we see 
that the 21.0 keV is only detectable at the l-cr level at 

^ |http: //xspec .gsf c .nasa .gov I 
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the start of the burst. The 14.0 keV is highly significant 
at the start of the burst, and remains detectable at the 
> l-cr level for ~130 s. The 8.8 keV feature is the only 
one that remains detectable at the > l-cr level beyond 
130 s from the burst onset. 

3.2. Pulse Profile Changes 

Many AXP outbursts a re accompanied by significant 
pulse profile changes (e.g. iKaspi et al.ir2003l ). To search 
for these in 4U 0142+61, for each observation, we gen- 
er ated 64-bin pulse profiles using the method described 
in iDib et all (|2007[ ). We then aligned the profiles with 
a high signal-to-noise template using a cross-correlation 
procedure in the Fourier domain. Then, for each group of 
observations in the active phase, we summed the aligned 
profiles, extracted the DC component from the summed 
profile, and scaled the resulting profile so that the value 
of the highest bin is unity and the lowest bin is zero. 

The resulting pulse profiles are shown in Figure [5] in 
chronological o rder . The panels in the left-most column 
are from Dib et al.l (|2007t ) and show pulse profiles in the 
6 years prior to the entry into the active phase. Notice 
the slow change in the height of the dip between the 
peaks. The panels in columns 2 and 3, marked with the 
corresponding letter in the top right corner, show pulse 
profiles for each of the data groups in the active phase 
that were marked with a letter in Figure [TJ The 4 plots 
in bold correspond to the observations marked with "x" , 
"bl", "b2", and "b3" in Figure [H 

The pulse profile evolution can be described as fol- 
lows: prior to the active phase, features in the double 
peaked pulse profile were evolving on a timescale of sev- 
eral years (see panels C6 to CIO corresponding to RXTE 
Cycles 6 to 10). Then, in the first observation of the ac- 
tive phase (panel x), the pulse profile became suddenly 
triple-peaked. It was also triple-peaked in the following 
observation (panel bl), in which a burst occurred. It re- 
mained multi-peaked for the following two groups of ob- 
servations (panels G and H). Then, an observation with 
multiple bursts occurred (panel b2). In that observation, 
also having a triple-peaked pulse profile, the middle peak 
was taller than the other two. Following this observation, 
the pulse profile seemed to be slowly recovering back to 
its double-peaked long-term shape (panels I, J, K, and 
L). Another burst observation interrupted this evolution 
seven months later (panel b3). In that observation, a 
large burst was detected. In the pulse profile of that ob- 
servation, the left-most peak was significantly taller than 
usual. The event that caused this change had apparently 
no effect on the following observation which occurred 2 
days later: the profile went back to being double-peaked 
(panels M and N). To summarize, the pulse profile be- 
came multi-peaked at the beginning of the active phase. 
Following the second observation with bursts, the profile 
started to recover to its double-peaked shape. The evo- 
lution was only temporarily interrupted for the duration 
of the third observation with bursts. The behavior of the 
pulse profile in the 2—4 and 4—10 keV bands was similar. 

Note that from Figure [5] alone we can compare the 
sizes of the peaks to each other, but we cannot track the 
evolution of the heights of each peak separately. In order 
to do that, we must scale the pulse profile of each group 
of observations by the average pulsed flux of that group. 
This analysis is presented in Section r3. 4. 31 



We also performed an analysis of the Fourier compo- 
nents of the pulse profiles. The results are shown in Fig- 
ure [S] The variations in the power of the fundamental 
Fourier component are shown in panel 1 , and that of the 
second harmonic in panel 2. Note how the amplitude of 
the fundamental varied monotonically prior to, but not 
during, the active phase. Also note how the power in the 
second harmonic was already back to its pre-active-phase 
level before the last burst occurred. 

3.3. Burst Rotational Phases 

An important factor in understanding the origin of the 
bursts is the rotational phase at which they occur. The 
phases of the bursts are shown in Figure [71 For each 
burst observation we created at 31.25-ms time resolution 
lightcurve and folded it using our timing solution (see 
§ 13. Sp . We then phase-aligned these folded profiles by 
cross-correlating them with the long-term pulse profile 
template. Our phase-aligned folded profiles are shown 
in Figure [7] (histograms in the last row). Superposed 
on each folded profile are two curves. The top curve is 
made of the 5 Fourier harmonics that best fit the his- 
togram. The bottom curve is made of the best-fit 5 har- 
monics after the removal of the 4 seconds centered on 
each burst. Note how the two curves in a given plot 
are similar, demonstrating that the additional peaks in 
the profiles are not due to the burst. The first three 
bursts occurred near the middle "new" peak in the pro- 
file. Burst 5 occurred near the old tall peak of the profile 
(see Fig. [5]). The phase of burst 6 corresponds again to 
a new peak in the profile, this time located where the 
previous small peak used to be. The coincidences of sev- 
eral bursts with new, transient profile features that are 
present even when the actual burst data are removed 
are suggestive of lower-level transient emission from the 
same physical location, with the burst being the extreme 
of this emission's luminosity distribution. 

3.4. Pulsed Flux Analysis 
3.4.1. Short-Term Pulsed Flux Analysis 

Previous AXP bursts are often accompanied by short- 
term pulsed flux enhancements (e.g. Gavriil et al. 2006). 
To search for these, for each of the three observations 
containing bursts, we made two barycentered time series 
in count rate per PCU, one for the 2—4 keV band and 
the other for 4—20 keV. We included only photons de- 
tected by PCUs that were on for the entire duration of 
the observation. The time resolution was 31.25 ms. We 
removed the 4 s centered on each burst from each time 
series. Then, we broke each time series into segments 
of length ^^500 s. For each segment, we calculated the 
pulsed flux using two different methods. 

First, we calculated the RMS pulsed fiux using 



RMS 



\ fe=i 



K.^+^b.^)), (2) 



where ak is the k^^ even Fourier component defined as 
Qk — jj^iLiPi '^'^^ i'^'^^^ / ^ '^ak is the uncertainty of 
Ofc, bk is the k^^ odd Fourier component defined as bk 
= jj'^^iPi sm{2TT ki/N), is the uncertainty ofbk, i 
refers to the phase bin, N is the total number of phase 
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TABLE 2 

4U 0142+61 Burst Temporal and Spectral Properties 





April 2006 




June 2006 




February 2007 


Parameter^ 


Burst 1 


Burst 2 


Burst 3 


Burst 4 


Burst 5 


Burst 6 


Burst day (MJD UTC) 
Burst start time (UT) 
Rise Time'' 

Phase'' 
Tqo Fluence" 
Tqo Fluencef 
Peak FluxS 
Peak Fluxli 


53831 
07:09.55.544(7) 
7+4 
'-3 

s.otlJ 

0.503(7) 

65±6 

16±3 
299±59 

75±19 


53911 
01:15:54.555(11) 
11+3 

0.487(4) 
1847±4 
229±11 
367dz65 
46±8 


53911 
01:15:55.119(43) 
40+22 

0.54(2) 
125±9 
26±3 
250±54 
53±12 

"'■0^-0.2?) 


53911 
01:16:9.216(4) 

4+3 
11 62+°-l5 
0.173(4) 
147±8 
38±3 
228±52 
58±14 


53911 
01:20:0.131(3) 
3+1 
92.540™5 
0.749(4) 
342±9 
43±4 
245±54 
31±7 

9 cr9+0.21 
^•0^_n.l3 


54138 
10:04:43.264(27) 
271^ 
3472.9^^4 
0.358(4) 
14377±41 
1455±53 
551±97 

253±22 
r, 19+0.08 

^■-^^-0.04 



^ All quoted errors represent l-cr uncertainties. 
^ Rise time (ms). 

" T90 duration (s), defined as the time it took to collect 90% of the total burst fluence. 

Phase (cycles), phases relative to the template shown in Figure [7] 

2-60 keV fluence in the T90 duration in units of counts PCU~^. 
' 2—60 keV fluence in the T90 duration in units of lO^^^erg cm"^. 
^ 2-60 keV peak flux for 64 ms in units of counts s~^ PCU~^. 
^ 2—60 keV peak flux for 64 ms in units of 10"'" erg s~^ cm"^ . 

' Blackbody temperature (keV). For burst 6, the temperature is that of the blackbody component after accounting for the apparent 
spectral features (see i|3. 1.21 1. 
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Fig. 2. — LEFT - 2—60 keV burst lightcurves binned with 31.25-ms time resolution as observed by RXTE. The red curves are the best-flt 
exponential rise and decay models (Eq. [TJ. RIGHT - 2-60 keV integrated burst lightcurves. For each burst lightcurve we subtracted the 
instrumental background, the background due to other sources in the FOV, and the contribution of the tails of neighboring bursts before 
integrating. Notice that some bursts have very long tails. The red curve is the integral of the model in the left panel. The blue curve is 
a step function whose height corresponds to the location of the total fluence as determined by extrapolating the best-fit integrated burst 
model (red curve) to infinity. 
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Fig. 3. — Panel lA: The spectrum of burst 6 extracted over 
a 5 s interval starting from 1 s before the burst peak. The red 
solid lines represents the best-fit model- a blackbody (blue curve) 
plus 3 Gaussian emission lines (green curves), all photoelectrically 
absorbed. The fit had a reduced of 1.2 for 32 degrees of freedom 
(DoF). Panel IB: Residuals, expressed in terms of the change in x^i 
after subtracting the best-fit model (blackbody -|- lines) shown in 
Panel lA. Panels 2A and 2B : Same, except here the best-fit model 
(red curve) consists of a blackbody (blue curve) plus 2 Gaussian 
emission lines (green curves), all photoelectrically absorbed. The 
fit had a reduced of 1-67 for 35 DoF. Panels 3A and 3B: Same, 
except here the best-fit model (red curve) consists of a blackbody 
(blue curve) plus a single Gaussian emission lines (green curve), 
both photoelectrically absorbed. The fit had a reduced of 1-89 
for 38 DoF. Panel 4A and 4B: Same, except here the best-fit model 
consists solely of a photoelectrically blackbody (red curve). The 
fit had a reduced of 2.54 for 41 DoF. 

bins, Pi is the count rate in the i^^ phase bin of the 
pulse profile, and n is the maximum number of Fourier 
harmonics to be tak en into account; here n=5 for consis- 
tency with lDib et~a l. (2007) and Gonzal ez etld] ([2009). 

While least sensitive to noise compared to other pulsed 
flux measurement methods, the RMS method returns 
a pulsed flux number that is affected by pulse profile 
changes. To confirm our pulsed flux results, we also used 
an area-based estimator to calculate the pulsed flux: 



1 



N 



-^Arca — ^ ^ Pmin ; 

i=l 



(3) 



where Pmin is the average count rate in the off-pulse phase 
bins of the profile, determined by cross-correlating with 
a high signal-to-noise template, and calculated in the 
Fourier domain after truncating the Fourier series to 5 



TABLE 3 
Burst 6 Spectral Fit 



Parameter 


Value'' 


Reduced x^(F)egrees of Freedom) 
2-60 koV Flux'' (10-3 erg s-'cm-^) 
2-20 keV Flux'' (lO"'' erg s-'cm-^) 


1.21(32) 

o qq+0.48 

0.»0_Q Jg 

2 67+0-12 
-0.06 


Blackbody Component 



Temperature, kT (keV) 
Radius'^ (km) 

2-60 koV Flux'' (10-3 erg s-'cm'^) 
2-20 keV Flux'' (lO"'-* erg s-'cm-^) 



1.8^ 



2.3 



+0.4 
0.3 

+ 1.0^/ d 
0.8 



iX 

0.53 
0.52 



— u.o 

(2.5 kpc) 

+0.06 
-0.02 
+0.02 
0.01 



Gaussian Emission Line 1 


Energy (keV) 
Width, 0- (keV) 

2-60 keV Flux'' (10"^ erg s-'cm-^) 
2-20 kcV Flux'' (10-3 erg s-'cm-^) 


0+0.2 

s.»_o 2 

9+0-3 
250+0-03' 
250+0-011 


Gaussian Emission Line 2 


Energy (keV) 
Width, cr (keV) 

2-60 keV Flux'' (lO"'-* erg s-'cm-2) 
2-20 kcV Flux'' (10-3 erg s-'cm-2) 


14 0+0-4 
9 Q+0.4 
^•'-'-0.4 
1 c;7+0.19 
l-°'-0.06 
1 56+0-07 


Gaussian Emission Line 3 


Energy (keV) 
Width, cr (keV) 

2-60 keV Flux'' (10-3 erg s^'cm^^) 
2-20 keV Flux'' (10-3 erg s-^cm-^) 


21.8+0-^ 
3 0+1-'^ 

1 filci+0-48 
i.D10_Q j^g 

360+0-016 

U.JDU_Q QQg 



All errors -were extracted at the 1 cr level, 
l' All quoted fluxes are unabsorbed. In th e fit, the column density 
was held fixed at Nh=0.64:XW'^^ cm~'^ l|Durant &: van KerkwiikI 
12009 ). 

The blackbody radius was calculated assuming a distance to the 
source of 2.5 kpc (HuUeman et al. 2004). 

harmonics. The results are shown for i^RMS in Figure [8] 
for i^RMS {Faich gives consistent results). Note the sig- 
nificant increase in the 4—20 keV pulsed fiux in the 2006 
June observation following the cluster of bursts. This 
increase is not present in 2—4 keV. Also note the signifi- 
cant rise and subsequent decay of the pulsed flux follow- 
ing the large 2007 February burst. The pulsed flux was 
sufficiently enhanced in these two observations that one 
can see individual pulsations by eye in Figure [7] in the 
two panels marked with a star, containing the raw burst 
lightcurves with 31.25-ms time resolution. 

3.4.2. Long-Term Pulsed Flux Analysis 

For each of the 94 analyzed observations, we created a 
pulse profile (in units of count rate per PCU) using the 
same procedure as in Section 13.21 Then we calculated 
the pulsed flux for each observation using Equations [2] 
and[3l Data from PCU were omitted because the long- 
term trend in the pulsed counts is not the same as that 
in the remaining PCUs, presumably due to the loss of its 
propane layer. Data from PCU 1 were omitted after the 
loss of its propane layer as well, on MJD 54094. 

We extracted pulsed fluxes in the 2—4 and 4—20 keV 
bands using both the RMS and area pulsed fiux method 
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Fig. 4. — Spectral evolution of burst 6. Top: Burst lightcurve binned with 1-s time resolution and extracted over the same energy 
band used for spectral fitting, 2-30 keV. Middle: Surface plot of the difference between the burst spectrum and a simple continuum model 
(blackbody) in terms of change in . Each vertical slice in the surface plot corresponds to a spectrum that was calculated over a 5-s long 
window. This window was translated across the burst in 1-s long steps. For each spectrum, a best-fit blackbody model was subtracted. The 
color wedge on the right maps the colors to the resulting change in . The contours indicate the regions where the spectrum differed by 
the model by l-tr (68.3% confidence, blue contour), 2-u (95.4% condfidence, green contour) and 3-(t (99.7%, yellow contour), respectively. 
Notice how the feature at 21.8 keV (see Table|3]l is detected at the l-c level only at the start of the burst. The feature at 14.0 keV is highly 
significant at the start of the bursts, and remains detectable at the >l-(7 level for ~130 s. Notice, that beyond ~130 s the only feature 
listed in Table [3] that remains detectable at the >l-cr level is the feature at 8.8 keV. Bottom: The insets are examples of individual 5-s 
long spectra used to generate the surface plot. These spectra were extracted at 0, 50 and 140 s from the burst onset, respectively. The red 
curves indicate the corresponding best-fit blackbody model. 



because of the numerous pulse profile changes around the 
times of the bursts. There is no evidence of a long-term 
change in the pulsed flux associated with the bursts ex- 
cept, possibly, for a hint of an increase in the 2—4 keV 
band roughly midway between the second and third ob- 
servations containing bursts. We also note that the 
4—20 keV pulsed flux of two of the observations con- 
taining bursts are significantly larger than the long-term 
average. Removing the bursts from these observations 
does not change this result. 

We performed the same analysis for individual obser- 
vations in 2—10 keV for an extended period of time. This 
is shown in panels 2 a nd 3 of Figure El This plot is an 
update to Figure 10 of lDib et al.r(|2007f) . 

In order to verify that trends seen in panels 2 and 3 
of Figure [9] are not an artifact of the evolution of the 
response of the detector, we calculated the pulsed flux 
in erg s~^ cm~^ using a method that took the evolution 
of the respon se int o account. The method is described 
in IDib et al.l (|2007f ) and takes into account spectral fits 



obt ained from ima ging data (in this case XMM data, 
see [Gonzalez et aLl l20091 to convert counts to energy for 
each combined set of observations. For this analysis, we 
used data from all PCUs; however, data from PCU 1 
were excluded after the loss its propane layer. Data from 
PCU were included because the response matrices used 
took into account the loss of its propane layer. 

The results are shown in panel 1 of Figure [9] The 
first 6 points, c o rrespo nding to RXTE Cycles 6—10, are 
from iDib et al.l (|2007( ). The second to last point is ob- 
tained by combining all observations that occurred be- 
tween bursts 1 and 2, but omitting observations contain- 
ing bursts. The observations that we included took place 
during the exponential recovery of the possible glitch (see 
Section [375)1 . The last point in panel 1 was obtained by 
combining the observations that occurred after burst 2, 
but again omitting those containing bursts. The obser- 
vations we included took place after the end of the expo- 
nential recovery of the possible glitch (sec Section 13. 5p . 

In the first of the two data points in the active phase, 
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Fig. 5. — Pulse profile evolution of 4U 0142+61. Left-most col umn: normalized average 2—10 keV pulse profiles in the 6 years prior to 
the entry into the active phase in chronological order (panels from IDib et aLll2007f) . Middle and right-most columns: normalized, averaged 
2—10 keV pulse profiles for each of the data groups marked by a letter in Figure [l] after the entry into the active phase. The 4 plots with 
bold labels correspond to the observations marked with "x" , "bl" , "b2" , and "b3" in Figure [T] 



the pulsed flux in erg cm^^ is 18±8% larger than 
in the pre-active p hase. This is co n sisten t with the in- 
crease reported in iGonzalez et al.l (|2009[ ) in the same 
energy range. A hint of this increase can be seen in 
panel 2 although it appears less significant. This discrep- 
ancy can be accounted for by t he fact that the spectrum 
changed: IGonzalez et al.l (|2009f ) reported a temporary in- 
crease in the spectral hardness in an XMM observation 
of 4U 0142-1-61 immediately following the bursts. 

3.4.3. Combined Pulse Morphology and Pulsed Flux Analysis 

In Section 13.21 we calculated the Fourier components 
of the aligned average pulse profiles. This gave us the 
relative amplitude of the pulse profile harmonics in each 
group of observations marked with a letter in Figure [TJ 
In Section 13.4.21 we calculated the pulsed flux for every 
observation. Here, we compute a weighted average of 
the pulsed flux for each group of observations using the 
flux points calculated in Section 13.4.21 We then recon- 
struct the profiles for each of the groups from the first 



five Fourier components (not including the DC), scale 
them by the average RMS pulsed fiux for that group, 
and add the necessary offset for the lowest point on each 
curve to be zero. This means that the resulting scaled 
profiles return the correct pulsed flux. The advantage of 
this analysis is that we can now trace the evolution of 
each of the peaks independently. 

The results are presented in Figure [TUl In panel (a), we 
show the average scaled pulse proflle per observation for 
observations in the year preceding the active phase. The 
profile is double-peaked. In panel (b), we show in bold 
the scaled profile for the observation containing burst 1 . 
We also show the scaled profiles in each of the groups 
of observations that followed the burst. The profiles are 
triple-peaked. In panel (c), we show in bold the scaled 
profile for the observation containing bursts 2, 3, 4, and 5. 
We also show the scaled profiles in each of the groups of 
observations that followed the bursts. Note how the in- 
crease in the pulsed flux in the observation containing the 
bursts is not only a consequence of the appearance of the 
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Fig. 6. — Top: evolution of the power in the fundamental Fourier component of the pulse profile of 4U 0142+61 in the 2—10 keV 
band. Bo t tom: Evolution of the power in the second harmonic. In both panels, the points with large horizontal error bars are from 
IDib et aLl II2007I ). The remaining points are obtained from groups of observations after the entry into the active phase. The three bold 
points correspond to the observations with bursts. The three dashed lines correspond to burst epochs. 

new peak, but a result of the increase in size of all three 
peaks. Also note the evolution of the pulse profile back 
to being double-peaked. In panel (d), we show in bold 
the scaled profile for the observation containing burst 6. 
We also show the scaled profiles in each of the groups 
of observations that followed the burst. Note again how 
the pulsed fiux increase is due to both peaks increasing 
in size. 

3.5. Timing Analysis 

Many AXP outbursts and active phases a re accompa 
nied by interesting timing anomalies (e.g. |Kaspi et aT 
20031 IWoods et all 120041 : [Israel et all [2007al : iDib et al' 

Next 



with 64 phase bins, were cross-correlated in the Fourier 
domain with a high signal-to-noise template created by 
adding phase- aligned profiles from all observations. The 
cross-correlation returned an average pulse time of ar- 
rival (TOA) for each observation corresponding to a fixed 
pulse phase. The pulse phase (jj at any time t can be ex- 
pressed as a Taylor expansion, 



<j>{t) = Mto) + Mt - to) + - tof 



6 



(4) 



20091) . Next we consider the timing behavior of 
4U 0142+61, which, as we show, also exhibits interesting 
evolution at the start of the active phase. 

For all our RXTE observations of 4U 0142+61, pho- 
ton arrival times at each epoch were adjusted to the 
solar system bary center using the position obtained by 
iPatel et al.l (|2003D from Chandra data. They were then 
binned with 31.25-ms time resolution. In the timing anal- 
ysis presented below, we included only the events in the 
energy range 2.5—9 keV, to maximize the signal-to-noise 
ratio of the pulse. 

Each barycentric binned time series was epoch-folded 
using an ephemeris determined iteratively by maintain- 
ing phase coherence; see below. Resulting pulse profiles. 



where v = 1/P is the pulse frequency, v = dv/dt, etc., 
and subscript '0" denotes a parameter evaluated at the 
reference epoch t = to. To obtain ephemerides, the TOAs 
were fitted to the above polynomial using the pulsar tim- 
ing software package TEMPO^. 

As explained above, our first set of TOAs was obtained 
by aligning the folded observations with a template pro- 
file using a cross-correlation procedure. In order to de- 
termine to what extent this set of TOAs was affected by 
the pulse profile changes that took place in the active 
phase, we generated two additional sets of TOAs (sets 2 
and 3). 

^ See jhttp://www. atnf.csiro.au/research/puls ar/tempo| 
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Fig. 7. — Phase analysis of the 6 detected bursts. Each column corresponds to an observation in which bursts were detected. Top row: 
200-s segments of the time series containing the bursts in the 2—20 keV band. The time resolution is 1 s. Middle row: 4-s segments of 
the time series containing the bursts. The time resolution 31.25 ms. Bottom row: aligned folds of the burst observations shown below 
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top curve is made of the 5 Fourier harmonics that best fit the histogram. The bottom curve is made of the best-fit 5 harmonics after the 
removal of the 4 seconds centered on each burst. The arrows indicate the phases of each burst. 



For set 2, we made the assumption that the location of 
the ofF-pulse, determined by finding the minimum in the 
pulse profile after smoothing, is not affected by the pulse 
profile changes. We then generated TOAs by aligning 
the off-pulse regions of the folded observations with the 
off-pulse region of the template and extracting phase dif- 
ferences. The last set of TOAs (Set 3) was obtained after 
aligning the tallest peak of each smoothed pulse profile 
with that of the long-term template. 

Timing residuals for all three sets of TOAs are shown 
in Figure [TT] using two different ephemerides. The linear 
ephemeris used in the top panel is that obtained by fitting 
the first set of TOAs in the regions indicated by the hori- 
zontal arrows, with a phase jump in between. The linear 
ephemeris used in the bottom panel is the ephemeris that 
best fit the first set of TOAs prior to the active phase. 
The same long-term trend is present in the first two sets 
of timing residuals, with more scatter in the second set. 
The difference in phase between each bold circle and the 
corresponding empty circle represents our uncertainty in 
determining a fiducial point on the pulsar. In the third 
set of residuals, the outlier points represent the observa- 
tions where the largest peak was no longer the right-most 
peak. Apart from the outliers, the same trend seen in the 
other two sets of residuals is present in the third set. 

From here on, we assume that the presence of the same 
trend in our three sets of residuals is an indication that 
the TOAs in the first set were not significantly affected 
by the pulse profile changes. We therefore have used the 



first set of TOAs in the remainder of this Section. 

Using the above assumption, the results of the tim- 
ing analysis are shown in Figure 1121 The results of a 
segmented timing analysis are in panels 1 and 2. The 
results of a long-term timing analysis are in panels 3, 4, 
5, and 6 of the same Figure. 

For the segmented timing analysis, we used the 70 
TOAs that had the smallest uncertainties and omitted 
the TOAs of the 3 burst observations as well as the TOA 
of the first observation in the active phase. We then 
divided the data into segments of similar pulse profiles. 
For each data segment we found a linear ephemeris us- 
ing TEMPO^. We plotted the resulting ephemerides in 
panel 1 of Figure [T^] with the uncertainties returned by 
TEMPO. The timing residuals are shown in panel 2. No- 
tice how the slope between bursts 1 and 2 is more neg- 
ative than the long-term average. 

We then fit a linear trend through all the pre-burst 
observations, and another linear ephemeris through the 
post-burst observations after the pulse profile had started 
to return to the double-peaked shape. We did not include 
the TOA corresponding to the last observation contain- 
ing bursts. The results are plotted in panel 3 of Figure [T2l 
with uncertainties, with the residuals in panel 4. Notice 
the difference in the slope in the two regions. In partic- 
ular, extrapolating the two ephemerides to a point be- 
tween their times of validity makes it seem like a sudden 

* For the segmented timing analysis we ran TEMPO in mode 1. 
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Fig. 8. — RMS pulsed flux within the observations containing bursts. Each column corresponds to one observation. In each column we 
show, descending vertically, the 1-s resolution lightcurve with the bursts indicated, the 2-4 keV RMS pulsed flux, and the 4—20 keV RMS 
pulsed flux. The dotted line in each of the pulsed flux plots shows the average of the pulsed fluxes obtained after segmenting and analyzing 
the time series of the observation immediately prior to the one shown. 



spin down, i.e. an "anti-glitch" occurred. 

Finally, we included all the TOAs and did one global 
fit. In order to provide a good fit to the TOAs at the 
onset of the active phase, we had to assume that a glitch 
occurred on MJD 53809, with the glitch model consisting 
of a permanent change in ly and j> and a frequency change 
Ai^d that recovered exponentially on a timescale Td, i.e., 

^ iyo{t) + Au + Ai/de-^*-*'''/^'' + AO {t - tg), (5) 

where voit) is the frequency evolution pre-glitch given by 
voit) ~ I'ito) + ^{t — to), Av is a instantaneous unrecov- 
ered frequency jump, Av^ is the frequency increase that 
decays exponentially on a time scale Td, tg is the glitch 
epoch determined by setting the phase jump to zero, and 
Av is the post-glitch change in the long-term frequency 
derivative. The values of the fit parameters are listed in 
Table m 

From Table [31 we can see that the total sudden change 
in frequency which happened at the onset of the active 
phase {Avtot — Av + Avd) is positive, and Avd decayed 
exponentially. At the end of the decay, the net effect 
of the recovered glitch was a negative Av. That is, the 
data suggest that after Vd had decayed, there remained 
net spin- down relative to the undisturbed ephemeris from 
before the active phase. 

The fit parameter Avd, from which we conclude a sud- 
den spin-up, is controlled primarily by the first few TOAs 
in the active phase. It is therefore possible that it is af- 
fected by pulse profile changes. However, as shown in 
Figure [m the TOAs in the active phase were not signif- 
icantly contaminated by pulse profile changes. Also, the 
segmented analysis shown in Figure [T^ clearly indicates 
an initial spin-up. 

In contrast to Avd, the fit parameter Av, from which 
we conclude a net spin down, is primarily controlled by 



TABLE 4 

Spin and Glitch Parameters for 4U 0142+61'' 



Parameter 


Value 


MJD start 


53481.268 


MJD end 


54235.319 


TOAs 


93'' 


v (Hz) 


0.1150920955(12) 


V (10-1"' Hz s-l) 


-2.661(9) 


Epoch (MJD) 


53809.185840 


Glitch Epoch (MJD) 


53809.185840 


Ai^ (Hz) 


-1.27(17)xl0-* 


A!> (Hz s-i) 


-3.1(1.2)xlO-i'5 


Ai/d (Hz) 


2.0(4)xl0-^ 


td (days) 


17.0(1.7) 


RMS residual (phase) 


0.0168 



^ Numbers in parentheses are TEMPO-reported Icr uncertainties. 

A single TOA was omitted due to the very poor signal to noise 
ratio in the corresponding observation. 



the TOAs outside of the active phase, which were cer- 
tainly not affected by pulse profile changes (Figure [T2l 
panel 3). 

We therefore conclude that the pulsar likely suffered a 
spin-up glitch near or at the start of the radiatively active 
phase, but that the glitch 'over-recovered' such that long- 
term, its net effect is a spin- down of the pulsar. 

4. DISCUSSION 

In this paper, we have described the timing, pulse pro- 
file, and pulsed flux behavior of 4U 0142-1-61, during its 
2006-2007 active phase, the first such episode yet stud- 
ied from this source. Specifically we have shown that in 
addition to a sudden departure from a slow, systematic 
evolution of the source's pulsed flux and pulse proflle, 
this AXP also suffered a signiflcant timing event that is 
best described as a sudden spin-up glitch, followed by 
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Fig. 9. — Long-term pulsed flux time series in 2 
for 4U 0142-f61. Panel 1: area pulsed flux in erg s~^ cm 
for combined obse rvations. The so lid line marks the 29±8% in- 
crease reported in lDib et al.l II2007IV Panel 2: area pulsed flux in 
counts/s/PCU for individual observations. Panel 3: RMS pulsed 
flux in counts/s/PCU for individual observations. All panels: the 
dashed lines mark the burst epochs. The points marked with stars 
are the pulsed fluxes of the observations containing bursts. 

a large decay of the frequency jump such that the net 
effect was a slow-down with A1//1' ~ —8 x 10^^. In- 
terestingly, the pulsed X-ray flux showed no significant 
change apart from immediately following b ursts, in con 
trast to other AXP radiative outbursts (e.g 
20031 iTam et al.l 120081: flbrahim et afl 12004 



Kaspi et al 



Muno et al 



2007( ). Further, we have reported on six SGR-like bursts 
from the source that occurred during this period, one of 
which was notable for its unusual spectrum, which was 
poorly fit by an continuum model. 

4.1. The SGR-Like Bursts from 4U 0142+61 

Five of the six bursts reported here for 4U 0142+61 
had fast-rise-slow-decay profiles, with tails much longer 
than the rise times. These morphologies are similar 
to the class of AX P bursts labeled as "Type B" by 
IWoods et all ()2005D . Type B AXP bursts have also been 
characterized by arriving preferentially at pulse maxi- 
mum, and are seen more o f ten in AXPs compa red with 
SGRs (jWoods et all [2005t iCxavriil et all [200I . They 
were suggested as being due to a sudden rearrangement 
of magnetic fiel d lines anchored in the crus t following a 
crustal fracture (jThompson fc Duncanlll995l ). as opposed 
to reconnection events in the upper magnetosphere (e.g. 
lLvutikovll20M) . argued as more likely for the shorter, 
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Fig. 10. — Average pulse profile per observation in the 2—10 keV 
band in several groups of observations reconstructed from the first 
five Fourier components, and scaled to return the appropriate 
pulsed flux, (a) average scaled pulse profile for observations in 
RXTE Cycle 10, the last RXTE Cycle before the active phase. 
(b) Bold curve: scaled pulse profile of the observation containing 
burst 1. Thin curves: average scaled pulse profile for each of the 
groups of observations that followed burst 1 (groups G and H in 
Figure [TJ. (c) Bold curve: scaled pulse profile of the observation 
containing bursts 2, 3, 4, and 5. Thin curves: average scaled pulse 
profile for each of the groups of observations that followed the 
bursts (groups I, J, K, and L in Figure[T]l. (d) Bold curve: scaled 
pulse profile of the observation containing burst 6. Thin curves: 
average scaled pulse profile for each of the groups of observations 
that followed burst 6 (groups M, and N from Figure HI. 
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Fig. 11. — Top: Timing residuals for all three sets of TOAs. Bold circles: timing residuals obtained using the first set of TOAs obtained 
by cross correlation in the Fourier domain. Empty circles: residuals obtained using the second set of TOAs obtained by aligning the 
off-pulse regions of the pulse profiles. Triangles: residuals obtained using the third set of TOAs obtained by aligning the tallest peak of 
the pulse profiles. The linear ephemeris used to produce all three sets of residuals is that which best fit the first set of TOAs in the regions 
indicated with the horizontal arrows (with an arbitrary phase jump in between). Bottom: Timing residuals for the same three sets of 
TOAs. The linear ephemeris used is that which best fit the first set of TOAs in the pre-active phase region. Both panels: the dotted line 
indicates the epoch of the first observation of the active phase. The dashed lines indicate the burst epochs. 



symmetric "Type A" bursts that show no preference for 
pulse maximum. However, it is interesting that in the 
4U 0142+61 bursts, in spite of their B-typc morphol- 
ogy, no clear preference for arrival at or near pulse max- 
imum was seen (see Table [2|)- This may blur somewhat 
the distinction between the putative types, although the 
morphological distinction remains clear. 

It is also notable that the vast majority of bursts seen 
from three of the four AXPs from which bursts have yet 
been detected, IE 1048.1-5937, XTE J1810-197and now 
4U 0142+61, had long tails following fast rises. The ex- 
ception todate is IE 2259+586 for which, during its 2002 
outburst, the minority (roughly one dozen out of 80) of 
the bursts seen were of this form, the majority being of 
Type A, sim ilar to those c lassic ally seen in SGRs, and 
suggested bv lWoods et al.l (|2005f ) to be magnetospheric. 
A possible hint regarding the origin of the different types 
may lie in that IE 2259+586 was mid-outburst when its 
bursts were observed, whereas the other sources' bursts 
all occurred in the days/months following the commence- 
ment of an active period, presumably as the pulsar re- 
covered from a major event, rather than mid-event. 



4.1.1. Spectral Features 

The complicated spectrum of burst 6 is puzzling. The 
most significant feature at ~14 keV is particularly in- 
teresting. Emission features at similar energies were ob- 
served from two out of the three bursts from IE 1048.1- 
5937 (jGavriiret al.ll2002l : IGavriil et al.H2006l) and in one 
out o f the four bursts from XTE J1810-197 ()Woods et all 
|2005| ). We have reanalyzed all these burst spectra con- 
sistently as for 4U 0142+61. In Figure [13] we plot the 
spectra of all AXP bursts with likely emission lines in 
their spectra. Notice that all have broad features that 
occur between 13 and 14 keV. There is possible evidence 
for features at --8 and --21 keV in the 4U 0142+61 burst; 
there may be hints of these features in some of the other 
burst spectra as well. 

If the features are interpreted as of a cyclotron origin, 
it is puzzling that three sources with spin-down-inferred 
magnetic fields that span a range of a factor of three pro- 
duce features at similar energies. The similarity in energy 
suggests a mechanism that is not very sensitive to the ex- 
act value of the magnetic field. Moreover, it is unclear 
why such features are seen only in some bursts. Thus far, 
all AXP burst spectral features have been observed with 
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Fig. 12. — Panel 1: Frequency vs. time plot of local linear 
ephemerides with uncertainties. The linear trend corresponding to 
the local ephemeris directly before the active phase is subtracted 
from all local ephemerides. Panel 2: Timing residuals after sub- 
tracting the ephemerides shown in panel 1. Panel 3: Frequency vs. 
time plot (with uncertainties) of the best-fit linear ephemerides in 
the pre-active phase region and in the post burst 2 region. The 
same trend as in panel 1 was subtracted. Panel 4: Timing resid- 
uals after subtracting the ephemerides shown in panel 3. Panel 5: 
Frequency vs. time plot (with uncertainties) of the best-fit long- 
term ephemeris which includes a glitch at MJD 53809 followed by 
a fast exponential recovery. The same trend as in panel 1 was sub- 
tracted. Panel 6: Timing residuals after subtracting the ephemeris 
shown in panel 5. 
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Fig. 13. — Burst spectra of all observed AXP bursts with signif- 
icant emission lines, as seen by RXTE. Panel lA: Burst spectrum 
of 4U 0142-(-61 burst 6. The dotted line indicates the continuum 
(blackbody) component of the best-fit model. Panel IB: Residuals, 
expressed in terms of the change in i after subtracting only the 
continuum component of the best-fit model. Panels 2 A and 2B : 
Same but for burst 4 of XTE J1810-197 Csee lWoods et al.ll2005l ). 
Panels 3A and 3B : Same, but for burst 3 of IE 1048.1-5937 (see 
IGavriil et al|[200^ . Panels 4A and 4B : Same, but for burst 1 of 
IE 1048.1-5937 (see IGavriil et al|[2003) . 

RXTE. Their absence from many other i?XT£^observed 
bursts, together with the absence of any response feature 
in the PCA near 13-14 keV, is supporting evidence for 
their veracity. Nevertheless, it is important to confirm 
these features with other instruments. Observations of 
AXPs during active phases using the Swift XRT or the 
large-area LAXPC aboard the planned Astrosat mission 
would be particularly useful in this regard. 

4.2. The Net Spin-Down Event in 4U 0142+61 

The timing glitch reported on in i i3.5l had recovery 
fraction, defined as Q = /\vd/{^Vd + Ai/), equal to 
1.07 ±0.02. Q > 1 implies that the net frequency change 
after the transients have decayed is negative (see panel 3 
of Figure [Hand Table [J). Indeed, the net Ai^/i^ ob- 
served is —1.1 ± 0.1 X 10"''. Such a Q > 1 has not 
previously been reported in any AXP but has been seen 
following magnetar-like r adiative behavior in one high-B 
rotat ion-powered pulsar ([Living stone et al.l[2009( ). Note 
that iMorii et al.l ([2005f ) reported a candidate glitch in 
1999 from this source. If such a glitch occured, based 
on the ephemerides reported in iDib et all ([2007[ ). the 
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fractional frequency change, after any recovery, would 
have been +1.9(2) x 10~'^-8.6(2) x 10"^ depending on the 
glitch epoch, but manifestly positive. 

In the standard model for glitches in rotation-powered 
pulsars, the neutron-star crust contains superfluid neu- 
trons that rotate faster than the bulk of the surrounding 
matter (see, e.g.,lAnderson fc Itohlll975l: [Anderson et al.l 



fTMllAlpar et al.lll984l:IAlpar fc Pine£lll993D . This is ar- 
gued to be a result of the fact that the external magnetic 
torque acts on the crust and coupled core components 
only, with the uncoupled superfluid unaffected. The su- 
perfluid's angular momentum resides in quantized vortex 
lines whose density is proportional to angular velocity. 
The vortex lines are suggested to be pinned to crustal 
nuclei, and suffer strong forces due to the angular ve- 
locity lag that builds between crust and superfluid. A 
glitch in this picture is a sudden unpinning and outward 
motion of vortex lines, with a transfer of angular mo- 
mentum from superfluid to crust. In a magnetar, unpin- 
ning may occur due to the strong internal magnetic field 
as it deforms the crust plastically or cracks it violently 
(|Thompson fc Duncani ri996l. 

In the 4U 0142+61 event we describe in §3.5( if the 
standard glitch model is roughly correct and applies here, 
then some regions of the stellar superfluid were origi- 
nally spinning slower than the crust. Then the tran- 
sient increase in frequency would be a result of transfer 
of angular momentum first from a more rapidly rotat- 
ing region, with a subsequent angular momentum drain 
from the c rust to a more s l owly rotating region. As 
argued by [Thompson et al.l (|2000D . regions of slower- 
rotating superfluid can occur in magnetars, because the 
superfluid vortex motions are governed not be spin-down- 
related forces, but by advection across the stellar surface 
by the deforming crust. Those authors show that the 
number of vortex lines per unit surface area of crust 
can increase or decrease depending on the crustal mo- 
tion relative to the stellar rotation axis. They invoke 
this possibility to explain a possible "anti-glitch" seen in 
SGR 1900+14. However, that event was orders of magni- 
tude larger than w hat we have observed in 4U 0142+61, 
with Aiy/iy 10~^ (|Thompson et al.]|2000r ). and also was 
likely associated wi th the extremely en ergetic flare seen 
on 1998 August 27 (|Mazets et al.lll999l ). 

Independent evidence for slow crustal deformations in 
4U 0142+61 that could result in the regions of slower- 
rotating superfluid required to explain the spin-down 
may come from the long-term flux and pulse p rofile evo- 
lution previously reported for this source bv iDib et all 
(|2007f ) and shown again here in Figures 5 and 6. Al- 
though such variations could also be magnetospheric in 
natur e, possibly due to twi s ting of the magneti c field 
(e..g. [Thompson et al.l l2002t iBeloborodov fc ThompsonI 
|2007[ ). pr oblems wit h this interpretation exist, as dis- 
cussed bv IDib et al.l (|2007f ). For example, if magneto- 
spheric, it is puzzling that most if not all of the changes 
in pulse profile are seen below 4 keV, with none above 
6 keV. On the other hand, slow evolution of the sur- 
face emission (e.g. lOzel fc Giiverll2007[ ). hence structure, 
is perhaps more consistent with low-energy evolution. It 
remains to be seen if detailed modeling can reproduce the 
sort of slow pulse profile changes we have observed (Fig. 
5), quite apart from the sudden appearance of new peaks. 



hence greater harmonic structure, near bursts. We note 
that the latter has also recentl y been seen durin g periods 
of activity in IE 1048.1-5937 (l^i!^ et al.l l2008D. 

Previouslv lKaspi et al.l ()2000f ) noted that the existence 
of AXP glitches having properties comparable to those 
seen rotation-powered pulsars was not in itself evidence 
for AXPs being magnetars. This is because the stan- 
dard glitch model applies regardless of the mechanism 
by which the crust slows down; in principle a glitch can 
occur in an accretion-powered neutron star, since in this 
case too there should develop a crust/superfluid angu- 
lar velocity lag. However only in the context of magne- 
tars has the poss ibility of more slowly ro tating superfluid 
been suggested ([Thompson et al.l |2000[ ) . The net spin- 
down in 4U 0142+61 may thus add to the already large 
amount of evidence against any accretion-powered origin 
for 4U 0142+61. 

4.3. The Active Phase of 4U 0142+61 and Other AXP 

Outbursts 

In many ways, the 2006-2007 active phase of 
4U 0142+61 is similar to other phases of activity seen 
in AXPs: it was punctuated by short bursts, pulse pro- 
file changes were seen, and it was accompanied by a sig- 
nificant rotational anomaly. However the 4U 0142+61 
activity is unique in one interesting way: the pulsar 
did not suffer a large, long-lived pulse flux increase at 
the beginning of the phase. Indeed its pulsed flux (Fig. 
9), apart from very near bursts, has remained relatively 
stable. We note that observations with focusing X-ray 
telescopes may reveal some ph ase-averaged flux varia- 
tions (see iGonzalez et al.l l2009t) ; indeed pulsed fraction 
has been shown to be inversely correlated with tot al flux 
in IE 1048.1-5937 (|Tiengo et al.ll2005l : pR!^^reraIl l20081 . 
rendering pulsed flux variations smaller relative to to- 
tal variations. Still, for 4U 0142+61, such a correlation 
would have to conspire to render the pulsed flux steady. 
This seems unlikely, however only focusing telescope ob- 
servations can rule this out. 

Other AXP radiative outbursts have been suggested 
as arising from large magnetospheric twists, with as- 
sociated magnetospheric currents returning to the stel- 
lar surface and heating it, resulting in increased X-ray 
emission from the source (e.g. Thomp son et al. [20021 : 



IBeloborodov fc Thomps"M3 l2007l : iBeloborodovl l2009( ). 
Such twists are thought to represent a release of mag- 
netic energy and helicity from the internally wound-up 
magnetic field. 

On the assumption that the pulsed flux is a reasonable 
proxy for the total flux for 4U 0142+61, here, we find no 
clear evidence for significantly increased X-ray emission 
on time scales longer than a few minutes. Thus any large- 
scale magnetospheri c twist scenario is probl ematic. As 
originally showed bv [Thompson et al] ([2OO20 . the X-ray 
luminosity due to returning currents reheating the sur- 
face in a significantly twisted magnetosphere in general 
are comparable to that produced from internal processes; 
this is clearly not observed in the 2006-7 active period of 
4U 0142+61. 

Instead, for this source, long-term evolution of the 
crust, driven presumably internally by field decay, and re- 
sulting in multiple unstable configurations though with- 
out any large scale magnetospheric twists, could result 
in sudden cracking and local rearrangements. This could 
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be accompanied by bursts and profile changes, as well as 
with vortex line shifting. Why such surface motion does 
not produce significant magnetospheric twists is, how- 
ever, puzzling, given that field lines are thought to be 
anchored in the crust; perhaps only for large motions do 
field lines become sufficiently twisted for enhanced X-ray 
e mission to be pro duced. 

iDib et all (|2008[ ) showed that other AXP glitches have 
been unaccompanied by radiative changes. For exam- 
ple, AXP IE 1841-045 has glitched multiple times yet its 
pulsed fiux remains very steady. Such radiatively "silent" 
glitches may be a result of internal activity that does 
not result in significant twisting of outer magnetosphere 
field lines, where as radiatively "loud" glitches could be 
those for which the magnetosphere is impacted. We note 
in fact that thus far, the data are consistent with all 
AXP radiative outbursts being accompanied by timing 
anomalies. This may also be true of SGRs, though in 
those cases, timing anomalies are harder to establish be- 
cause of the difficulty in achieving phase-coherent timing 
in quiescence. 

5. SUMMARY AND CONCLUSIONS 

We have reported on anomalous X-ray pulsar 
4U 0142+61 entering an active phase which was preceded 
with by a long-term increase in pulsed fiux. The active 
phase, which commenced in 2006 March, consisted of a 
timing anomaly that can be described as a net "anti- 
glitch," that is, a net spin-down following an initial spin- 
up that decayed on a time scale of 17 days. Following 
the glitch, we detected six bursts from the pulsar, the 
first ever observed from this source. Despite 10 yr of 
RXTE monitoring, the bursts all occurred in the nar- 
row time span between 2006 April 6 and 2007 Febru- 



ary 7. The sixth and largest burst had a unusual spec- 
trum that cannot be described by any simple continuum 
model. Rather, it can be well fit if a prominent spec- 
tral feature at ~14 keV is included, as has been reported 
in other AXP bursts, as well as weaker features at '^8 
and ^^21 keV. The pulse profile of the source changed 
from double- to triple- peaked near the bursts, and un- 
derwent considerable evolution otherwise. At the burst 
epochs the relative intensity of the three peaks signifi- 
cantly varied. The pulse profile is now relaxing to its pre- 
active phase morphology. Most aspects of 4U 0142-|-61's 
emission changed during the active phase, with the no- 
table exception of the pulsed flux (except near bursts). 
This argues against this event being associated with a 
sudden magnetospheric twist, as has been invoked for 
other AXP activity, and is suggestive of crustal evolu- 
tion driven internally by the large magnetic field, though 
without significant magnetospheric twisting. We suggest 
that other, radiatively silent AXP glitches have a simi- 
lar origin, whereas radiatively loud AXP timing events 
occur when the crustal motions cannot avoid significant 
twisting of the magnetic field lines. 
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